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Summary 


The  investigations  on  this  program  have  included  the  measurement 
of  multiwavelength  scintillation  statistics  over  a  very  long,  highly  turbulent 
path;  the  analytical  and  experimental  study  of  finite -transmitter  and  wander¬ 
tracking  effects  on  target  ir radiance;  and  the  study  of  the  detailed  nature  of 
microthermal  turbulence  fluctuations,  including  turbulent  intermittency,  as 
related  to  short-term  scintillation  statistics.  The  latter  two  areas  represent 
ongoing  efforts,  and  are  reviewed  in  this  report. 

Finite -beam  experiments  are  being  conducted  with  a  wander -cancelling, 
tracking -transmitter  over  a  1,  6  km,  uniform  path.  A  new  prediction  is 
given  for  the  transverse  amplitude  correlation  length  under  near  field  con¬ 
ditions  in  strong  turbulence.  A  useful,  ^nysically  derived  approximation  to 
combined  diffraction,  wander,  and  short-term  beamspread  angles  or  mean 
target  irradiance  is  presented,  and  shown  to  agree  well  with  more  rigorous 
numerical  calculations;  experiments  under  strong  turbulence  conditions 
show  a  substantial  and  surprising  advantage  for  wander -tracking.  The 
normalized  variance  of  the  fluctuating  irradiance  (strength  of  fading)  due  to 
scintillation,  wander,  and  coherent  fading  is  also  discussed,  and  probability 
distributions  and  spectra  of  wander  and  irradiance  are  presented. 

Microthermal  fluctuations  due  to  turbulence  are  considered,  and  it 
is  shown  experimentally  that  the  statistics  of  appropriately  averaged  first- 
increments  or  derivatives  behave  in  accordance  with  recent  theoretical  pre¬ 
dictions.  Short-term  turbulence  and  scintillation  statistics  are  discussed, 
including  intermittency  effects.  It  is  pointed  out  that  the  1 .  vel-crossing  and 
conditional  probability  problems  for  the  scintillating  irradiance  are  compli¬ 
cated  by  the  non-Markovian  nature  of  this  process. 
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I.  Introduction 

Work  on  this  program  has  been  concerned  with  three  remaining 
aspects  of  turbulence -induced  scintillations  which  have  been  poorly  under¬ 
stood: 

1)  Scintillation  effects  far  into  the  multiple  scattering  or 
"saturation"  regime. 

2)  The  effects  of  beam -wave  or  finite -aperture  transmitters, 
including  the  influence  of  combined  wander,  scintillation, 
and  short-term  spread  on  target  irradiance.  This 
includes  the  advantages  and  parameter -dependencies  of 
wander -cancellation  through  tracking. 

3)  The  nature  and  effects  of  the  intermittency  of  turbulence 
on  scintillations. 

The  results  and  status  of  these  respective  efforts  will  be  reviewed 
in  the  following  sections.  The  latter  two  topics  represent  ongoing  investiga¬ 
tions  for  a  follow-on  program,  and  the  emphasis  in  the  present  report  will 
be  on  investigations  conducted  during  the  final  quarter  of  the  contract  period. 
The  reader  is  referred  to  preceding  Quarterly  Progress  Reports  for  details 
of  the  preceding  efforts. 

II.  Long -Path  Scintillations 

Scintillation  statistics  were  measured  at  4880  X  and  10.  6  jj,  wave¬ 
lengths  over  a  very  long,  low  path,  in  order  to  experience  very  high  inte¬ 
grated-path  turbulence.  The  pathlength  L  was  6  km.  The  results  are 
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presented  in  an  earlier  report  and  recent  paper;  these  include  the  observa¬ 
tion  of  saturation  of  scintillations  at  10.  6  |a,  and  heretofore  unobserved 
anomalies  in  the  behavior  of  covariance  functions  and  scintillation  spectra. 
We  briefly  discuss  here  some  recent  additions  to  the  understanding  of  this 
regime. 

3 

In  addition  to  Brown’s  numerical  analysis  for  a  two-dimensional 

case,  there  now  exist  two  highly  physical  treatments  of  multiple  scattering 
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or  scintillation  in  strong  turbulence.  '  These  treatments  are  generally 
similar,  although  the  results  disagree  as  to  the  behavior  of  the  log  ampli¬ 
tude  variance  with  increasing  pathlength  or  turbulence  "beyond  saturation!1.  ^ 
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Significantly,  both  treatments  predict  the  profound  change  in 
covariance  behavior  which  was  observed  in  the  long-path  experiments.  l>  2 
As  the  first-order  or  Rytov  theoretical  log  amplitude  variance  (<r  2) 
increases,  the  covariance  is  predicted  to  drop  rapidly  and  develop  a  long 
tail,  such  that  cwo  (i.  e.,  small  and  large)  transverse  amplitude  correlation 
scales  emerge  in  place  of  the  familiar  (  XL)1/2.  The  large  correlation 
scale  implies  poor  receiver  aperture  averaging  of  scintillations,  while  the 
short  scale  implies  the  existence  of  surprisingly  high  scintillation  frequen¬ 
cies  (several  kHz  at  visible  wavelengths);  both  of  these  phenomena  are 
observed. 

The  latter  effect  is  of  especial  importance,  since  scintillation 
experimenters  have  long  assumed  that  one  kilohertz  represents  an  ade¬ 
quate  frequency  response  for  measurement  apparatus,  and  that  the  neces¬ 
sary  condition  for  "point"  performance  of  a  small  r  reiver  was  that  its 
diameter  be  smaller  than  (XL)1/2.  R  is  now  clear  that  these  are  not  suf¬ 
ficient  conditions,  and  that  further  experiments  are  desirable. 

Although  our  10.6  p  results  were  not  affected,  the  interesting 
regime  of  »  1  should  be  further  investigated  at  shorter  wavelengths- 

For  example,  the  behavior  of  experimental  log  amplitude  variance  (  n  2  , 

"  't  cast  in  doubt;  although  we  suspect  that  the  negative  slope  or 
"supersaturation"  is  a  fact,  the  actual  slope  and  existence  of  an  asymptote 
are  not  verified. 

Yura  has  extended  his  analysis  to  explicitly  predict  scintillation 
spectra  in  this  regime.  We  have  just  received  this  work  as  of  this  writing 
and  we  will  compare  these  predictions  with  our  data  below  1  kHz.  It  is  clear 
that  the  general  features  are  in  agreement.  We  are  also  considering  the 
possibility  of  correcting  those  <r ^  values  which  are  affected  by  spatial  and 
electronic  filtering.  However,  Yura's  results  are  for  a  plane  wave  source 
and  omit  certain  coefficients;  it  would  be  far  preferable  to  conduct  new 
exje-iments  over  our  long-path  facility  using  very  small  receivers  and 
large  measurement  bandwidths. 


III.  Finite -Beam  Effects 

As  described  in  the  preceding  report,  the  equipment  which  has 
been  fabricated  for  measuring  finite -beam  effects  (including  tracking - 
out  of  atmospherically  induced  wander)  has  been  installed  at  our  regular 
field  site,  which  is  flat  and  uniform  and  includes  micrometeorological 
instrumentation.  The  tracking  transmitter  is  mounted  on  a  massive  con¬ 
crete  slab  which  is  isolated  from  the  transmitter  shack.  Wander -tracking 
and  target-ir radiance  signals  are  recorded  for  digital  processing.  During 
the  recent  period,  certain  minor  improvements  were  completed  and  the 
facility  is  now  under  full  operation  for  a  program  of  quantitative  experi¬ 
ments  at  6328  X.  High-turbulence  results  are  being  obtained  during  the 
summer,  and  lower -turbulence  data --which  will  scale  to  the  10.  6  p. 
case--will  be  conducted  through  the  coming  months. 

Efforts  are  continuing  on  the  clarification  and  unification  of  the 
theory  using  the  approach  first  reported  in  Ref.  1.  In  the  present  report, 
further  theoretical  considerations  are  discussed,  and  preliminary,  repre¬ 
sentative  experimental  results  are  given. 

A.  Amplitude  Correlation  Scale  (Scintillation  Patch  Size) 

The  nominal  size  of  the  transverse  amplitude  correlation  length 
(pa)  has  been  a  matter  of  considerable  confusion.  We  believe  that  it  is  now 
possible  to  confidently  predict  this  parameter  under  a  variety  of  conditions, 
as  follows. 


We  confine  ourselves  to  a  statistically  homogeneous  (horizontal) 
path,  and  denote  the  following  parameters: 

Optical/infrared  wavenumber  k 

6  ? 

Turbulence  strength  C  ^ 


Pathlenpth 
Coherence  scale^ 
Transmitter  diameter 
Rytov  log  amplitude  variance 
Beamwidth  at  receiver 


L 

P 


-  (k2LC  V3/5 
o  n 


2.c 

1  n 


w  (>  p  ) 
—  a 
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In  all  near-field  (b> ^ ^  )  conditions,  the  transmitter  is  assumed  to  be 
focused.  We  omit  coefficients  of  order  unity,  and  hence  we  do  not  dis¬ 
tinguish  between  p^  for  plane  and  spherical  wave  sources. 

Considering  scintillations  only,  i.  e.  omitting  beam  wander  effects, 
the  amplitude  correlation  scale  values  are  given  in  Table  I.  F^r  comparison, 
the  nominal  beam  sizes  are  also  given  (Ref.  1). 


Table  I 

Amplitude  Correlation  Scale  (p  )  and  Beamwidth  (w)  at  Receiver  Plane 

a 


Discussion 


Condition  (a)  corresponds  to  a  far -field  source  with  first-order 
theory  applicable.  ° 

Case  (b)  is  for  a  near-field  source  with  first-order  theory  applicable, 
and  is  the  only  condition  for  which  the  patch  and  total  beam  sizes  are  nominally 
the  same.  This  represents  the  "transmitter  smoothing  of  scintillations"  pre¬ 
dicted  by  the  first-order  theory  ari  illustrated  by  the  photograph  of  Fig.  1  of 

Ref.  9:  there  is  a  relatively  stable  central  spot,  with  little  scintillation 
around  it. 

Condition  (c)  represents  the  interesting  case  in  which  the  first-order 
theory  breaks  down  for  a  finite  transmitter  (b  >  Pq)  even  though  it  is  valid 
for  a  plane  or  spherical  wave  source  <  1).  This  is  the  realm  predicted 
in  Ref.  9  and  illustrated  by  the  photograph  of  Fig.  3  in  that  paper,  which 
shows  the  proliferation  of  diffraction-scale  patches.  It  has  also  been  pre¬ 
dicted  by  Goshelashvily.  ^ 

Conditions  (d-f)  are  not  describable  by  first-order  theory.  There 
is  every  physical  and  dimensional  reason  to  believe  that  the  results  of  Yura's 
analysis  for  a  plane  wave  (Ref.  4),  i.  e.  that  ^  becomes  Pq  for  <r  2  >  l, 
will  also  apply  to  a  far -field  source.  The  result  given  fo/conditJn  (f),  for 
a  focused,  large  source,  is  derived  from  simple  asymptotic  (or  breakpoint) 
consistency  with  cases  (c)  and  (e).  This  is  a  novel  result,  and  says  that  the 
effect  of  a  large,  focused  transmitter  aperture  (b)  in  a  strong,  multiple - 
scattering  regime  is  to  yield  a  scintillation  patch  size  =  p  X  (p  /b)(«  p  ). 
This  physically  corresponds  to  further  breakup  of  ^-size/patches  into 
extremely  small  spots,  it  is  to  be  hoped  that  Yura's  analysis  can  be  extended 
to  cover  such  a  condition,  which  realistically  applies  for  a  typical  near -field 
(focused)  condition  with  visible  wavelengths  and  strong  turbulence. 

A  particularly  interesting  behavior  may  be  illustrated  as  follows. 

Suppose  that  we  have  a  near-field  transmitter  and  wish  to  determine  the 
behavior  of  pa  as  turbulence  increases.  We  will  then  progress  from  case  (b) 
tc  (c)  and  then  (f),  as  shown  in  TabJe  IA: 


Table  IA 


Amplitude  Correlation  Scale  (p  )  for  Near-Field  Transmitter  and  Increasing 

3i 

Turbulence 


Condition 


^  b  < 
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kb 
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kb 
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This  predicts  that  the  patch  size  will  remain  that  of  the  transmitter  diffrac- 

l/2  2 

tion  scale  until  p^  decreases  to  (L/k)  7  (i.  e.  tr  =1),  and  then  it  will  decrease 

as  p  ^/b  as  turbulence  strength  increases. 

In  all  cases,  the  scintillation  spectra  will  have  a  bandwidth  ~  v  /'p  , 

n  sl 

where  is  the  transverse  wind  velocity.  Case  (f)  therefore  implies  very 
high  scintillation  frequencies  (>  10  kHz)- -enhanced  by  a  factor  of  b/p  compared 
to  those  observed  in  Ref.  1  and  predicted  in  Ref.  7. 

We  remark  here  that  this  predicted  phenomenon  of  an  even  greater 
degree  of  breakup,  for  strong  and  increasing  turbulence  and  a  focused,  near¬ 
field  transmitter,  may  be  manifested  by  an  increase  in  log  amplitude  variance 
or  total  scintillation  level  over  that  naively  expected;  that  is,  the  final  asymp¬ 
totic  behavior  in  Figure  20  in  Ref.  1,  which  is  the  log  amplitude  variance  for 
large  b,  may  be  larger  or  more  complicated  than  expected.  Of  course,  the 
high  frequencies  and  small  patch  sizes  involved  make  accurate  scintillation 
measurements  difficult  and  in  many  cases  may  indicate  diminishing  practical 
importance  in  target  illumination  situations. 

Finally,  we  note  that  the  ultimate  breakdown  of  turbulence  MTF 


theory  (and  e.  g.  the  existing  expressions  for  p  )  for  stroig  turbulence 

scattering  is  not  known.  Nevertheless,  the  present  expressions  for  p 

Ka 

represent  a  significant- -and  probably  sufficient- -extension  from  the  first- 
order  scintillation  theory. 

Experiments 

The  predictions  for  cases  (d-f)  in  Table  I  may  be  verified  through 
photographs  of  the  scintillation  patches  at  the  receiver  plane,  as  was  done 
for  cases  (b)  and  (c)  in  Ref.  9.  This  will  be  undertaken  if  sufficiently 
straightforward  with  existing  facilities.  Ultimately,  it  would  be  highly 
desirable  to  employ  large -bandwidth  electronics  with  very  small  receivers, 
to  accurately  measure  covariance  and  variance  effects  in  this  regime. 

Case  (d)  then  represents  the  situation  of  Sec.  II,  and  case  (f)  a  large -aperture 
extension  of  it.  Note,  however,  that  a  very  long  path  is  not  required  for 
case  (f). 

B.  Finite-Beam  and  Tracking  Effects  on  Mean  Target  Irradiance 

In  this  section,  we  extend  the  description  begun  in  Ref.  1  of  finite - 

transmitter  and  wander -tracking  effects  on  the  first  moment  of  target 

irradiance  (average  illumination)  on  the  long-term  beam  axis.  The  approach, 

which  utilizes  recent  developments  in  Huygens -Fresnel  theory11,  12  and 
•x  H-13  . 

reciprocity,  involves  the  simple  mean-square  addition  of  the  individual 

angular  beamspread  mechanisms,  and  resolves  contradictions  and  incon¬ 
sistencies  in  the  literature.  As  will  be  seen,  the  expressions  represent 
an  excellent  approximation  to  more  rigorous  numerical  results,  while  retaining 
clear  physical  meaning  and  dimensional  consistency.  The  discussion  will  also 
point  out  theoretical  work  which  remains  to  be  completed.  The  present  dis¬ 
cussion  will  be  confined  to  the  case  of  a  negligible  inner  scale  of  turbulence; 

the  extension  to  non  zero  f  is  indicated  in  Ref.  1. 

o 

We  first  point  out  that  integral  expressions  for  the  long-term  (no 
wander -tracking )  and  short-term  (wander-tracked)  mean  irradiance  are 
given  in  Ref.  14,  and  that  these  expressions  can  be  numerically  evaluated 
for  an  arbitrary  transmitter  wave.  An  exact  analytical  solution  is  probably 
not  possible,  which  indicates  that  the  approximate  solution  given  below  in 
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terms  of  separate  physical  mechanisms,  simply  fails  to  take  into  account 
some  degree  of  "coupling  1  between  these  mechanisms.  Numerical  results 
for  a  gaussian  beam  may  be  readily  obtained  and  would  be  q\ute  useful. 
Related  numerical  results  for  a  uniformly  illuminated  lens  have  been  given 
by  Ron,  and,  through  reciprocity,  by  Fried.  1^)'  17  The  latter  will  be 
discussed  further  below  and  related  to  the  present  approximate  analysis. 

1.  Calculation  of  Spread  Angle  from  Physical  Beamspread 


Mechanisms 

The  separate  beamspreading  or  irradiance -reducing  mechanisms 
are  readily  derived  to  within  a  numerical  constant.  We  initially  defer  the 
specification  of  these  constants,  which  avoids  the  specification  of  the  trans¬ 


mitter  beam  wave  (gaussian,  truncated  gaussian,  uniform  circular,  etc.), 
thus  preserving  generality.  Also,  it  avoids  the  question  of  competing 
angular  definitions  such  as  l/e  diameter,  second  spatial  moment,  etc., 
which  can  yield  significantly  different  results.  18'  19  It  is  implicit  in  our 
development  that  angular  spreads  are  defined  in  terms  of  on-axis  irradiance, 
which  is  the  important  quantity. 

We  now  write  the  expressions  for  each  respective  beam-spreading 
mechanism  as  follows: 

Diffraction  Angle 

Mean-Square  Diffraction  Angle  ~  (la) 

where  k  is  the  optical/infrared  wavenumber,  and  b  the 
transmitter  beam  diameter  by  an  appropriate  measure. 
Geometric  Angle 

Mean-Square  Geometric  Spread  Angle  ^  b  (  —  -  -  )  (lb) 

D  R 

where  L  is  the  transmitter-target  distance  and  R  is  the 
transmitted  wavefront  curvature.  In  particular,  this  term 
is  zero  for  the  focused  case  (R=L). 

Wander  Angle 

We  obtain  the  wander  angle  by  applying  the  reciprocity  principle  to 
receiver  image  dancing  (Ref.  1): 


9- 


Mean-Square  Wander  Angle  ~ 


D0(b)  ,  5/3 

2  2  "  (b/po)  ■■  1 . 

k  b  ,  2  2 

k  b 


(Po/b> 


1/3  1 


v2  2 
k  p 

ro 


where  is  the  phase  structure  function,  and  p 

is  the  coherence  scale,  as  observed  at  the  receiver, 

11-13 

for  a  conceptual  point  source  at  the  target. 

It  may  easily  be  shown  that  this  is  a  geometrical-optics  expression,  which 
does  not  actually  involve  k.  However,  the  final  form  shown  is  convenient 
for  comparison  with  other  terms. 

Atmospheric  Beam  Spread 

Finally,  we  consider  the  short-term  atmospheric  beam  spread. 

This  has  previously^'  ^  been  given  in  the  asymptotic  form  for  b  »p  , 

2  2  ° 
as  (l/k  p^  );  this  is  the  regime  of  "beam  break-up".  However,  the  details 

of  this  mechanism  for  b  comparable  to  p  can  be  important,  especially  if 

wander  is  cancelled  out,  and  we  write 

Mean-Square  (Short-Term)  Spread  Angle  ~  f(b/ p  )  — ~ — -  (Id) 

O  .  L  Ld 

k  p 

o 

where  f(b/p  )  can  be  obtained  from  a  complete  (numerical)  analysis  and  will 

be  discussed  later.  Note  that  f( —  -*■  ®)  =  1. 

po 

Total  Angular  Spiead 

We  now  assume  a  focused  condition,  and  write  the  total  mean-square 
2 

angle  <J>  from  Eqs.  (la-ld)  as 

2 

<3?  =  (diffraction)  +  (wander)  +  (short-term  spread) 


a  m 

kV  +  2 


/b_\ 
\  po< 


+  c3  Hf 


2  2 


In  the  event  that  wander  is  tracked  out,  C^=0. 

For  the  case  of  fixed  pQ  and  fixed  transmitter  size  b,  respectively, 
Eq.  (2)  can  be  usefully  written  as 
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2,  2  2  1 

*  k  po  >—  )2  + 
o 


-1/3 


2,  2  2 

$  k  b 


=  ci  + 


[*£ 


t  c  f- 

3  tpo; 


+  C3  ft 


Asymptotic  log  plots  of  these  expressions  are  given  in  Figs.  la,b.  In 
Fig.  la,  increased  abscissa  corresponds  to  increased  transmitter  size, 
while  in  Fig.  lb,  it  corresponds  to  decreasing  p  or  increasing  turbulence 
strength  and/or  pathlength.  The  short-term  spread  function  f(b/p  )  cannot 
usefully  be  broken  down  into  asymptotic  linear  segments,  and  is  shown 
qualitatively  in  Figs.  la,b. 

Finally,  with  the  determination  of  C  in  Eq.  (2),  the  complete, 
continuous  curves  of  $  will  be  specified.  The  conceptual  curves  with  and 
without  wander  tracking  are  shown  in  Figs.  la,b. 

2.  Angular  Parameter  a 

As  discussed  below  in  conjunction  with  the  second  moment  of 
irradiance,  a  parameter  of  importance  is  the  ratio  of  mean-square  wander  to 
total  short-term  spread.  From  Eq.  (2),  this  is  given  by 

c  P  "5/3  b"1/3 

a  =  —L _ ° _  (3) 

ci  r  1 

_ 1  C  f  —  — * 

o  +  3  (p  /  P  C 

b2  vo/  ^o 

For  a  small  aperture  (b«  p  ),  diffraction  predominates  and  we  have 
b5//3 

a  ^  TR  ~  D<b>  »  14a) 


where  D  is  the  wave  structure  function.  For  a  large  aperture  (b»p  ),  we 

.  ./P„V/3 


/P  \ 
Q  q  __ 2  I 

lb  / 
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3.  Corresponding  First  Moment  or  Mean  Irradiance 
The  mean  target  irradiance  at  the  long-term  beam  centroid 
is  simply  proportional  to  the  reciprocal  of  $2.  Hence, 

T  =  i 

o  2  _  2 
$  L 


where  Iq  is  a  measure  of  the  irradiance  at  the  transmitter.  For  reference 
in  later  discussion  below,  we  also  write  this  as 


4.  Comparison  with  Numerical  Results;  Determination  of  Con¬ 
stants 

For  a  given  transmitter  wavefront,  the  constants  in  Eq.  (5a) 
might  in  some  cases  be  determined  from  a  more  detailed  consideration  of  the 
physical  beamspread  mechanisms.  Alternatively,  they  may  be  determined 
through  comparison  with  complete  numerical  solutions.  Such  a  solution  is 
available  from  the  analysis  of  an  optical  heterodyne  receiver;17  through  the 
principle  of  reciprocity,  13  this  signal  has  a  one-to-one  correspondence  with 
the  target  irradiance  for  a  uniformly  illuminated,  focused,  transmitting 
aperture. 

The  heterodyne  analysis  includes  results  with  and  without  wave - 
front  tilt-tracking,  which  correspond  respectively  to  the  cases  of  wander - 
tracking  and  a  static  transmitter.  In  particular,  the  maximum  improvement 
in  mean  irradiance  to  be  obtained  through  proper  control  of  aperture  size  (b) 
and  the  employment  of  wander -tracking  (C2  =0)  can  be  predicted. 
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We  first  note  that,  in  Eq.  (5a),  we  have  one  redundant  constant. 
We  rewrite  the  expression  as 


1  po 


+  f/i 


where  I  is  the  asymptotic  target  irradiance  for  large  b,  and 
2  2  °° 

depends  on  k  Pq  (see  Eq.  2a).  We  then  note  that  the  intersection  of  the  dif¬ 
fraction  and  large-  b,  short-term  beam  spread  asymptotes  (Fig.  1)  occurs  in 
the  heterodyne  analysis  at  b  =  r  =  2.  0986  p  *.  Hence  we  have 


Ci  P. 


(2.0986  p  )* 

o 


or  Cj  =  4.  40 


for  a  uniformly  illuminated  transmitter. 

To  determine  the  detailed  function  f(b/pQ),  we  use  the  heterodyne 
results  for  the  tracking  case  (C^  =0).  Identifying  i/l^  with  "antenna  gain  G"  ; 
defined  in  Ref.  17,  we  write  from  the  reciprocal  of  Eq.  (6): 


I  4.40 
G  ’  •’"/Po2 


=  H 


The  curve  f (b/pQ)  as  obtained  from  the  numerical  results  of  Refs.  16  and  17 
is  given  in  Fig.  2.  For  each  value  of  the  abscissa,  the  function  f  may  be 
interpreted  as  the  coefficient  on  the  short-term,  atmospherically  induced 
beam  spread  relative  to  the  large-aperture  value  (l/k2?  2). 


*Note  that  this  number  was  incorrectly  given  as  2.  15  in  Ref.  1.  Also,  the 
appropriate  pQ  for  the  present  case  is  the  spherical  wave  value.  1 
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For  b/po  <  3,  the  physical  interpretation  of  this  function  is  less 
simple  since  it  increases  for  decreasing  b,  while  the  wavefront  distortion 
and  associated  spreading  (with  linear  tilt  removed)  should  approach  zero. 
This  occurs  because  f  in  Eq.  (8)  becomes  the  small  difference  between 
two  large  numbers;  in  fact,  this  fu  -ction  actually  manifests  the  "coupling" 

between  diffraction  by  b  and  P().  As  ±  further  decreases,  f  has  a 
negligible  effect  on  G.  ° 

For  large  b/pQ,  f  is  seen  to  approach  unity  surprisingly 
slowly.  This  regime,  will  be  discussed  further  in  Section  III-5, 

,  N°W  that  C1  and  f  have  been  determined,  we  may  choose 

C2  such  that  Eq.  (6)  will  best  approximate  the  numerical  results  for  the 
no-tracking  case.  An  approximate  optimum  value  is 

C  '=  1.80 

2  ’  (9) 

which  establishes  the  relative  weighting  of  the  beam-wander  effects  on 
mean  irradiance  for  this  transmitter  configuration. 

The  usefulness  of  the  present  approximation,  i.  e.  u  sing 
explicit  physical  mechanisms,  can  be  tested  by  comparing  the  curve  of 
Eq.  (6)  with  the  numerical  results  of  Refs.  16  and  17  for  the  no-tracking 
case.  This  comparison  is  shown  in  Figure  3,  where  the  tracking  case 

(C2  =0)  iS  als°  shown  for  completeness.  The  approximation  is  seen  to  be 
quite  good,  which  indicates  that  wander  is  largely  independent  or  uncoupled 
from  the  other  mechanisms.  The  relationship  between  these  curves  and  the 
(reciprocal)  asymptotic  curves  of  Fig.  la  is  readily  apparent. 

The  case  corresponding  to  Fig.  lb,  i.  e.  with  fixed  b  and 

variable  Pq,  can  also  be  plotted  using  the  same  constants  C  '  and  C.  '  and 

function  f(b/po).  The  ordinate  must  then  be  renormalized  by  multiplying 

the  denominators  of  Eq.  (6)  by  b2/pQ2  (see  Eq.  2b);  the  normalizing  quantity 

or  asymptote  at  small  b/P<)  then  corresponds  to  the  irradiance  for  weak 

turbulence,  which  is  determined  by  transmitter  (free -space)  diffraction  and 
depends  on  kb. 
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5.  Concept  of  Short-Term  Coherence  Scale  (p  )  and  Related 

0 

Approximations  st 

14 

In  a  recent  paper,  Yura  discusses  a  "short-term"  atmospheric 
muti  ax  coherence  function  (MCF)  in  which  linear  wavefront  tilt  (or  wander) 
has  been  removed.  16  The  expression  for  this  MCF  is 

r  5/3 r  v  l/3 

Mst(p)  =  00) 

where  p^  is  the  coherence  scale  used  in  the  preceding  expressions,  and  rep¬ 
resents  the  l/e  point  of  the  conventional  or  "long-term"  MCF.  Yura  then 
similarly  defines  a  "short-term"  coherence  scale,  p^  ,  as  the  l/e  point  of 
Mst(P);  ^  may  be  noted  that  p  is  a  function  of  the  (uniformly  illuminated) 
transmitter  aperture  size  b.  the  purpose  of  the  present  discussion  is  to 
show  that  certain  approximations  are  incorrectly  applied  in  Ref.  14,  and  that, 
contrary  to  that  development,  the  physical  utility  of  Pq  is  marginal  at  best. 
Hence  that  paper  does  not  properly  treat  the  problem  Stof  wander  cancella¬ 
tion  and  beam  spread. 

From  Eq.  (10),  we  determine  the  value  of  Po  ^  as  a  function  of 

st 

h/ p^  &s  follows.  F rom  the  definition  of  p  ,  we  have 

st 


The  solution  of  this  equation  is  plotted  in  Figure  4.  It  should  be  noted  that, 

for  (b/pQ)<  1.21,  Pq  does  not  exist;  i.  e. ,  M  does  not  have  a  1 /e  point. 

st  st 


In  Ref.  14,  the  author  approximates  the  solution  of  Eq  (11)  by 


st 


■1/3 


(12) 


This  is  also  plotted  in  Figure  4,  from  which  it  can  be  seen  that  the  approxi¬ 
mation  is  poor  for  b/p  <  10.  Unfortunately,  in  Ref.  14  the  approximation 

o 

is  utilized  extensively  for  1.  0  <  b/p^  <_  co,  and  in  particular  in  the  regime 

below  10,  (Note-  the  author's  abscissa  is  the  reciprocal  of  that  used  here). 

Therefore,  Figs.  2-4  in  that  reference  are  without  value  and,  in  fact,  p^ 

does  not  even  exist  for  p  /b>  0.  83. 

o 


st 


A  related  difficulty  with  Ref.  14  is  that  the  author  equates  the 

z  2  2 

short-term  (or  wander-removed)  atmospheric  beam  spread  to  1/k  p  , 

s  t 

which  is  equivalent  to  identifying  the  function  f  of  Eq.  (Id)  and  Fig.  2 
However,  there  is  no  physical  basis  for  this  use  of  p 


with  p  2/p  2 
o  r  o 

s  t 

except  for  b  »  p  ,  in  which  case  p  «  p  and  p 

st  st  st 


become  s 


st 


trivial.  Similarly,  he  equates  the  long-term  (including  wander)  beam  spread 
2  2 

to  l/k  p  ,  which  is  likewise  incorrect  except  for  b  »  p  .  Hence,  even  if 


the  author  had  utilized  the  correct  values  of  p  ,  the  interpretations  in  terms 

°st 

of  angular  spread  or  mean  irradiance  would  be  erroneous. 

In  order  to  gain  further  insight  into  this  problem,  we  may  define 


Pq  such  that  atmospherically  induced  beam  spread  is  given  by  l/k2p 


for 


eff 


any  value  of  b/p q.  From  Eq.  (2): 


eff 


1 


1 


,2  2 
k  p 

°eff 


v2  2 
k  p 

m 


(i) 


-1/3 


+  f  — 


(13) 


where  C  -  C^/C^  -  1*80  in  the  long  term  case  discussed  previously,  and 

C  =  3  in  the  short  term  or  wander-tracked  case.  In  the  long  term  case, 

P o  *s  aPProximatecl  by  p^  only  for  b/p^  large  enough  that  wander  is  negligible 

anc  f  converges  to  unity.  For  the  short  term  case,  Pq  is  approximated 

by  P  only  for  b/p  large  enough  that  f  converges  to  (p  ^/p  2  ). 
st  °  ro  "o 

Finally ,  we  note  from  Eq.  (12)  that  st 
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The  approach  to  unity  is  thus  very  slow,  due  to  the  -l/3  exponent;  it  is 
within  one  percent  for  b/ p^  >  4x10^.  Hence  the  convergences  in  the 
preceding  paragraph  occur  very  slowly, 

6.  Use  of  Approximate  MCF's 

The  basic  Huygens -Fresnel  integral  expression  for  mean  irra- 

14  20 

diance  has  been  solved  analytically  for  a  gaussian  beam  by  Titterton, 

using  an  approximate  MCF  in  which  the  5/3  exponent  on  (p/p  )  is  replaced 

by  2.  This  has  been  shown  to  be  equivalent  to  eliminating  all  wavefront 

21 

distortion  terms  except  linear  tilt.  Unfortunately,  this  seemingly  minor 
approximation  in  the  exponent  leads  to  results  which  are  dimensionally 
inconsistent  and  therefore  misleading.  However,  the  approach  is  worthy 
of  discussion. 

Specifically,  the  result  corresponding  to  Eq.  (5b)  is 


I  =  I 


,2,  2 

k  b 


1  +  2 

P  * 

ro 


An  earlier  analysis  of  wander 


alone 


22 


gives  the  result 


(15) 


(15a) 


I 


1 

1  +  2q 


(16) 


where  a  is  the  ratio  of  mean-square  wander  to  instantaneous  beam  width. 
By  comparing  Eqs.  (15)  and  (16),  it  is  tempting  to  identify  a  with  b2/p  2, 
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and  to  identify  the  two  denominator  terms  of  Eq.  (15)  as  diffraction  and 
wander,  respectively.  However,  wander  is  in  fact  a  geometrical-optics 
phenomenon,  and  comparison  with  Eq.  (5b)  shows  that  the  final  term  in 
Eq.  (15)  is  dimensionally  related  to  atmospherically  induced  beam  spread 
rather  than  wander.  The  dimensionally  correct  expressions  for  a  are 
given  in  Eqs,  (3,4). 

Similarly,  direct  calculations  of  wander  from  the  definition  of 

Ref.  23,  or  from  the  phase  structure  function  and  reciprocity  as  applied  to 
20 

wavefront  tilt,  lead  to  the  same,  incorrect  results  when  the  approximate 
MCF  is  employed. 

This  approximate  MCF  approach  will  be  discussed  further  in 
connection  with  the  second  moment  in  a  subsequent  section. 

C.  Finite-Beam  and  Tracking  Effects  on  Irradiance  Fading  or 
Second  Moment 

1  .  Analysis  of  Fading 

As  discussed  in  Ref.  1,  the  physics  of  fading  may  be  separated 
into  the  following  three  mechanisms: 

a)  Wander -fading  • 

b)  Scintillation  (first-order) 

c)  Beam-breakup  or  coherent  fading. 

The  statistics  of  wander  fading  are  predicted  to  be  beta-distributed,22’24 
while  the  latter  two  mechanisms  are  thought  to  be  log  normal.  To  determine 
the  total  fading,  the  individual  variances  are  added. 

Wander  Fading 

25  22 

Utilizing  results  from  Esposito,  Titterton  has  shown  that  the 
normalized  variance  of  irradiance  due  to  wander  is  given  by 


Z 

^wander 


4a  +  1 


(17) 


where  a  has  been  defined  previously  and  is  correctly  given  by  Eq.  (3). 


Scintillation 


First  order  scintillation  for  a  finite  transmitter  has  been  solved 

by  several  investigators.  ^  Detailed  numerical  results  for  the  log  amplitude 

Z  6 

variance  are  given  by  Kerr  and  Eiss.  In  the  present  context,  where  we 

are  unfortunately  forced  to  deal  with  the  linear  irradiance  rather  than  the 

27 

more  meaningful  log  quantities,  it  is  customary  to  write 


2 

scin  i’lation 


(18) 


In  this  expression,  tr  is  the  log  amplitude  variance  for  a  point  transmitter, 
and  prepresents  a  "transmitter  aperture  smoothing"  term  which  is  related 
to  the  reciprocal,  receiver-aperture  smoothing  in  a  heterodyne  receiver. 

In  order  to  write  Eq.  (18),  the  scintillations  must  be  assumed  to  be  log 

9 

normal.  Providing  that  the  transmitter  wavefront  is  carefully  adjusted,  the 

26 

value  of  p  can  decrease  significantly  as  the  aperture  is  increased.  Rather 

than  relate  p  to  the  approximate  analysis  of  direct-detection  receivers  in 

the  literature,  as  is  customary,  we  may  precisely  define  p  by  specifying 

2 

that  the  log  amplitude  variance  results  in  Ref.  26  are  equal  to  (  p  <r  ). 

Beam  Breakup  or  Coherent  Fading 

As  the  transmitter  aperture  b  becomes  greater  than  p  ,  the  beam 

o 

at  the  target  breaks  up,  i.  e. ,  it  shows  substantial  short-term  spread.  This 
is  accompanied  by  an  increase  in  fading  with  increased  aperture  size  which 
is  not  predicted  by  first-order  scintillation  theory,  and  which  is  related  to 

27 

the  reciprocal  case  of  "coherent  fading"  in  an  optical  heterodyne  receiver. 

The  effect  is  presumably  an  added  mechanism  for  log  normal  scintillations. 
Unfortunately,  this  situation  has  not  been  well  analyzsd;  relevant  work 
includes  the  approximate  coherent  fading  analysis  of  Ref.  27,  and  the  finite  - 
transmitter,  multiple -scattering  analysis  of  Ref.  10.  The  latter  purports 
to  directly  show  the  breakdown  of  first-order  theory  for  b  >  p  . 

'v  r  O 

For  an  arbitrarily  large  transmitter,  the  results  of  Ref.  27  imply 

that  fading  will  increase  without  limit.  However,  this  is  not  experimentally 
9 

observed.  Rather,  for  the  large,  focused-transmitter  case,  we  expect  an 


-19- 
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asymptote  on  the  order  of  the  plane  wave  variance.  As  discussed  above  in 
Sec.  Ill -A,  the  situation  may  be  especially  complicated  when  the  Rytov 


c  T»  1. 


28 


Brown  has  given  insight  into  the  large -aperture  problem  by 

assuming  a  set  of  N  mutually  incoherent  radiators  on  the  aperture 
2  2 

(N=  b  /p  »  1),  and  showing  that  the  normalized  variance  is 


CP. 


1  - 


large  aperture 
coherent  fading 


b2/p02 


(19) 


where  the  quantity  in  parentheses  is  the  irradiance  variance  for  a  spherical 

wave  source  (Eq.  18).  This  shows  that  the  asymptote  is  unity,  which  is 

2 

reached  at  smaller  b/p  when  cr  is  large. 

In  general,  we  must  admit  that  the  coherent  fading  regime  is  poorly 
understood,  and  we  therefore  represent  this  term  by 


g 


(20) 


coherent  fading 

where  g  is  a  function  whose  general--but  not  detailed — behavior  is  known. 
Combined  Fading 

If  we  now  assume  the  independence  o.-  ;he  fading  mechanisms,  we 
write  from  Eqs.  (17,  18,  and  20): 


2  4a2  43  <r2 

T  i.;  ,45T  +  (e  ’»  +«  • 

total 


(21) 


The  proper  independent  variable  for  a  is  (b/p  ),  as  can  be  seen  from  Eq.  (3) 
by  multiplying  the  denominator  and  numerator  by  b  .  The  first-order  scintil- 

2  _  9 

lation  term  depends  on  cr  and  b/ s /  JU/k  ,  and  we  note  that  cr  is  related  to  p 

o 

by 

—  5/3 


cr  =  0.  228 


N/T7k 


(22) 
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The  g -dependence  is  unknown  but  probably  depends  primarily  on  (b/p  ), 

2  _  ° 

and  secondarily  on  <r  and  perhaps  \l  L  /k.  Hence,  a  single,  universal 

curve  cannot  be  drawn. 

The  asymptotic  behavior  of  the  wander-fading  term  vs  (b/p  )  is 

obtained  from  Eqs.  (3)  and  (17),  and  is  discussed  in  Ref.  1.  For  the 

wander -tracked  case,  a  is  taken  as  zero. 

2 

The  general  behavior  of  cr^.  vs  b  is  shown  in  Fig.  5a  for 

2  total 

fixed  p  ,  \l  L/k,  and  c r .  To  derive  a  similar  curve  for  variable  o  and 
o  rQ 

fixed  b  and  n/ I./k  (and  hence  fixed  p),  we  use  Eq.  ( 22)  to  write  the  first- 
order  scintillation  term  as 


4(3  <r 

e  r  -  1  = 


exp  |  0.  912  (3 


[n/L Zk 


-"15/3 


-1 


«  0.912  (3 


QL2fl5/3 


(23) 


where  we  assume  (3  «  1.  The  resultant  curve  is  shown  in  Fig.  5b. 

2.  Comparison  with  Approximate  Analyses 
Using  the  Huygens -Fresnel  approach,  the  approximate  MCF 
discussed  in  a  previous  section,  and  certain  other  approximations,  Titterton 
has  derived  the  following  expression  for  the  normalized  variance:20 


(24) 


It  is  interesting  to  compare  this  result  with  Eq.  (21);  if  we  identify  a  with 
2  2  1 
(b  / pQ  ),  we  may  identify  the  first  term  with  wander-fading  and  the  second 

with  scintillation.  There  is  no  term  for  coherent  fading,  and  a  cross -coupling 
of  wander-fading  and  scintillation  has  apparently  arisen.  However,  as  dis¬ 
cussed  previously,  this  identification  of  a  is  clearly  inconsistent.  Some  of 
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the  same  approximations  were  also  employed  in  Ref.  27.  It  is  obviously 
very  desirable  to  obtain  a  (probably  numerical)  result  which  is  free  of  these 
approximations. 

D.  Representative  Experimental  Data 
1.  Introduction 

In  this  section  we  present  some  preliminary,  representative 
data  on  finite-beam  and  wander -tracking  effects.  These  data  were  taken 
under  high  turbulence  conditions,  such  that  b»  p^.  Extensive  additional 
data  will  be  collected  under  a  variety  of  conditions. 

A  summary  of  the  data  being  collected  appears  in  Table  II, 
and  the  experimental  and  data  processing  parameters  are  given  in  Table  III. 


Table  II.  Experimental  Data  Collected  for  Finite-Beam  and  Wander- 
Tracking  Effects 

Tracker  (transmitter)  end:  Real-time  wander  signal  (servo-drive  signal) 

Target  (receiver)  end:  rieal-time  linear  and  log  irradiance  signal 

Mean  irradiance 

2 

Strength  of  turbulence  (C  ,  p  ) 

n  o 

Real-time  microthermal  fluctuations 
Meteorological  parameters 


Table  III.  Experimental  and  Data  Processing  Parameters 
Transmitter  and  beacon  wavelengths:  6328  X 
Transmitter  power:  75  mW 
Path:  1.  6  km,  flat  farm  land 
Beam  height:  1.8  m 

Mean-irradiance  averaging  time:  100  sec 
Variance  of  log  and  linear  irradiance  and  wander  signal: 
Digital  computation  with  10  sec  averaging  time, 

300  sec  data  records 

Irradiance  and  wander  spectra:  Digital  computation  with 
0.  1  Hz  resolution,  100  sec  averaging  time;  analog 
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computation  with  1  Hz  resolution,  100  sec  averaging 
time 

Irradiance  probabilitv  distributions:  digital  computation 
Digital  sampling  rate:  1  kHz 
Target  receiver  bandwidth:  1  kHz 
Target  receiver  aperture:  0.  6  cm  ( 

Target  beacon  aperture:  0.  6  mm 
Transmitter  aperture:  truncated  gaussian  beam,  precisely 
focused,  15.2  cm 

Inner  scale  determination:  from  spectrum  of  microthermal 
signal 


The  tracker  servo  resonance  is  at  several  hundred  Hz,  and  the  tracker 
closed-loop  response,  which  is  a  function  of  gain  settings  and  atmospheric 
conditions,  is  much  faster  than  the  beam  wander  effect. 

The  tracker  system  utilizes  a  quadrant  photodetector  and 
(x*y)  galvanometric  mirror  scanners  to  maintain  the  target-beacon  image 
coincident  with  the  axis  of  the  transmitter  laser;  hence,  the  instantaneous 
mean  angle -of -arrival  of  the  beacon  wavefront  at  the  transmitter  aperture 
is  matched  by  the  outgoing  transmitter  wave.  Since  the  beacon-image --and 
transmitted  beam- -are  spread  out  by  the  atmosphere,  the  system  actually 
operates  on  the  centroid  of  the  image  distribution.  It  has  been  pointed  out29 
that  this  centroid  differs  from  the  conventional  first-moment  definition, 
owing  to  the 'nature  of  the  quadrant  detector: 

Tracker  "centroid"  displacement  =  J"i(x)  sgn  (x)  dx  (25) 

where  i(x)  is  the  image  irradiance  distribution.  However,  it  is  physically 
apparent  that  the  distinction  is  not  critical,  and  this  contention  is  supported 
by  the  detailed  experimental  results. 

A  much  more  important  consideration  is  the  following.  In 
our  application  of  reciprocity,  in  which  we  use  a  fast-tracking  transmitter  and 
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target  to  cancel  beam  wander,  it  is  implicit  that  the  entire  transmitted  beam 
is  refracted  uniformly  by  the  atmosphere:  that  is,  the  image  of  the  point- 
beacon  must  accurately  indicate  the  tracking  deflection  necessary  to  cancel 
the  wander  or  centroid  displacement  on  the  target.  Under  high  turbulence 
conditions  (b  »  Pq),  the  beam  is  substantially  broken  up,  and  rays  reaching 
the  transmitter  from  the  beacon  may  not  accurately  indicate  the  refraction 
of  more  than  part  of  the  beam. 

This  situation,  which  is  similar  to  an  isoplanatic  patch  argument, 
can  result  in  inaccurate  centering  of  the  instantaneous,  transmitted  beam 
centroid  on  the  target  beacon.  That  is,  the  instantaneous  centroid  of  the 
amorphous,  broken-up  beam  may  m^ve  around  rapidly;  in  fact,  the  beam 
can  then  be  envisioned  in  terms  of  separate  bundles  which  wander  independently. 

This  effect  may  be  expected  to  average-out  over  relatively 
short  periods  (e.g.  ,  a  fraction  of  a  second),  which  is  confirmed  by  visual 
observation  at  the  target.  Furthermore,  the  significant  wander,  even  for  a 
badly  broken-up  beam,  tends  to  be  surprisingly  uniform  over  the  beam, 
which  indicates  the  prism-like  effects  of  turbulence  scales  which  are  larger 
than  those  included  in  the  theoretical  (inertial  subrange)  model.  The  latter 
phenomenon  accounts  in  part  for  the  surprising  improvement  in  target 
illumination  which  is  observed  with  wander -tracking  under  conditions  of 
strong  turbulence,  as  discussed  in  a  later  section. 

Perfect  centering  of  the  transmitted  beam  centroid  on  the 
desired  target  point  can  in  principle  be  achieved  with  a  sophisticated,  image- 
processing  tracking  system  which  operates  on  the  extended  image  of  the 
finite,  illuminated  target  area. 

2.  Wander  Signal 
RMS  Wander  Angle 

The  recording  of  the  wander -tracking  signal  at  the  trans¬ 
mitter  was  begun  just  prior  to  the  writing  of  this  report.  In  the  first  such 
run,  the  results  for  the  rms  angular  wander  were  (perhaps  fortuitously)  in 
good  agreement  with  the  theoretical  prediction: 


Measured  C  2 :  5.  14X10"14m"2^3 
n 

Calculated  p  :  6.  5  mm 
ro 

Approximate  effective  aperture  (b)  =  15.2  cm 
Approximate  predicted  rms  wander  angle  = 


Theoretical  predictions  exist  for  the  power  spectral  behavior 
of  the  optical  phase  difference  between  two  points.  3°’  31' 32  This  may  algo 

be  applied  to  the  spectrum  of  image  dancing  or  wander.  32  For  the  particular 
data  run  under  discussion  here,  the  wind  velocity  was  very  low,  and  the 
applicable  prediction  is 

Power  spectrum  of  phase  difference  •>.  u‘8/3  .  (26) 

The  power  spectrum  and  frequency-weighted  power  spectrum 
of  the  wander -tracking  signal  is  shown  in  Fig.  6.  The  log-log  slopes  wore 

very  nearly  equal  to  -8/3  and  -5/3  respectively,  which  is  consistent  with 
the  above  prediction. 

Wander  Distribution 

It  is  generally  assumed  that  a  one -dimensional  component 
of  beam  wander  is  normally  distributed.  22  The  distribution  of  the  vertical 
component  of  the  wander  signal  for  the  data  run  under  discussion  here  is 
shown  in  Figure  7.  Marked  asymmetry  about  the  mean  is  noted,  which  is 
probably  a  manifestation  of  slow,  large-scale,  vertical  refractive  effects. 

3.  Target  Signal 

We  now  present  representative  results  for  the  target  irra- 
diance  with  and  without  wander -tracking,  from  eight  data  runs  in  strong 
turbulence.  These  runs  represent  the  condition  (b  »  pj,  and  in  some  cases 
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may  be  affected  by  the  inner  scale,  since  is  small.  Although  the  inner 
scale  has  been  determined  through  microthermal  spectra  for  each  run,  we 
have  not  yet  attempted  to  modify  the  theoretical  predictions  for  those  cases. 
Although  we  are  taking  advantage  of  summer  conditions  for  the  high- 
turbulence  data,  the  more  interesting  case  of  b  comparable  to  p  will  be 
examined  later;  in  particular,  that  is  the  case  which  scales  to  10  micron 
systems. 

Irradiance  vs  Time 

The  target  irradiance  on  the  beam  axis  is  shown  in  Fig.  8a 
without  wander-tracking  and  Fig.  8b  with  tracking.  The  irradiance  is 
smoothed  over  a  1  second  averaging  time,  in  order  to  reduce  the  influence 
of  scintillations.  From  these  records,  for  which  b/pQ  =  34.  9,  the  substan¬ 
tial  advantage  in  mean  irradiance  and  reduced  fading  to  be  obtained  from 
tracking  is  qualitatively  obvious. 

Mean  Irradiance 

The  mean  irradiance  on  the  bean,  axis,  corrected  for  changes 

in  transmitter  power  from  run  to  run,  is  shown  in  Fig.  9.  Since  b  is  fixed 

\ 

and  po  is  varied,  this  case  corresponds  to  the  reciprocal  of  Fig.  lb;  we 
expect  a  minus  two  (log-log)  slope  al  .arge  b/p  ,  although  Fig.  2  demonstrates 
that  the  slope  will  still  be  distorted  by  f(b/p  )  at  b/p  <  300.  For  variable 
Po*  we  do  not  as  yet  have  a  detailed  numerical  prediction  such  as  that  of 
Fig.  3  for  variable  b. 

The  ratio  of  the  mean  irradiance  for  tracking  vs  no  tracking  is 
shown  in  Fig.  10.  The  theoretical  prediction,  which  is  a  universal  curve,  is 
calculated  from  the  results  of  Ref.  16  for  a  uniform  transmitter  aperture. 

It  is  seen  that  the  improvement  realized  by  tracking  is  generally  significantly 
greater  than  that  predicted  fox1  such  large  b/ PQ i  we  surmise  that  this  is  due 
to  beam  refraction  from  turbulence  scales  which  are  larger  than  the  outer 
scale  and  hence  not  included  in  the  theory. 

Fading 

The  log  amplitude  variance  of  the  fluctuating  target  signal 
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is  given  for  the  wander -tracking  and  no-tracking  esses  in  Table  IV,  as  a 
function  of  measurement  bandwidth.  In  the  no-tracking  case,  the  total 
fading  consists  of  components  due  to  beam-breakup  scintillation  (or 
coherent  fading),  and  wander,  while  in  the  tracking  case,  the  wander  is 
removed.  Since  wander  is  a  low-frequency  phenomenon,  the  reduction 
in  variance  is  seen  to  be  most  substantial  when  higher  frequency  (scin¬ 
tillation)  components  are  either  spatially  or  electronically  filtered  out. 

The  full -bandwidth  (1  kHz)  data  of  Table  IV  may  be  related 

to  the  curves  of  Fig.  5b  by  determining  the  normalized,  linear  irradiance 
2 

variance  ( tr^.  )  instead  of  the  log  amplitude  variance.  The  latter  quantity 
is  more  indicative  of  dynamic  range  (the  dynamic  range  of  fading  in  dB  = 
100  cr  for  the  log  normal  case1),  but  the  limited  theoretical  treatments 
available  for  this  regime  tend  to  deal  with  c r*.  For  the  log  normal  case, 

which  does  not  apply  when  wander -fading  is  substantial,  the  two  variances 
are  related  by 


2  4cr 
°I  =C 


(27a) 


„-2  1  i  ,  2 

o-  =  4  In  (<rl  + 


(27b) 


Rather  than  making  use  of  these  relationships,  we  will  re¬ 
process  the  data  directly  fort t*.  It  may  be  noted  that,  according  to  Eq. 
(27a),  o-j  exceeds  unity  for  the  full -bandwidth  data  of  Table  IV. 

Spectra  of  Fading 

Typical  power  spectra  for  the  log-amplitude  with  and  without 
tracking  aie  shown  in  Fig.  11,  and  the  low-frequency  contribution  from 
wander  is  readily  apparent.  The  frequency-weighted  spectra  will  often 
show  a  double -peaked  behavior  for  the  no-tracking  case,  where  the  low 

and  high  frequency  peaks  correspond  respectively  to  wander  and  scintilla¬ 
tions  . 

Probability  Distribution  of  Fading 

The  probability  distribution  of  fading  without  tracking  includes 
the  combination  of  log  normal  scintillations  and,  theoretically,  exponentially 


Table  IV.  Log  Amplitude  Fading  Variance  vs  Bandwidth,  With  (T)  and  Without  (NT)  Wander -Tracking 


distributed  wane or-faaing.  These  distributions  have  not  yet  been  measured 
in  detail;  the  two  fading  mechanisms  may  be  largely  separated  through 
spectral  filtering,  as  discussed  above. 

Typical  distributions  of  log  amplitude  are  shown  in  Figure  13. 

As  expected,  the  distribution  becomes  much  more  closely  log  normal  when 
wander -tracking  is  employed. 

IV.  Short-Term  Turbulence  and  Scintillation  Statistics;  Intermittency 

A.  Introduction 

In  this  effort  we  are  concerned  with  modeling  the  microthermal  field 
in  such  a  manner  that  short-term  scintillation  statistics  can  be  predicted 
and  related  to  turbulence  statistics.  The  goal  in  this  work  is  to  understand 
measurement  confidence  intervals,  conditional  scintillation  statistics, 
scintillation  level  crossings  and  deep  fading,  and  ultimately  phase  and 
imaging  effects,  as  related  to  finite  or  short  measurement  times  and  to  the 
often  observed  sporadic  (intermittent)  nature  of  the  turbulence  mechanism. 

In  order  to  implement  this  study,  we  are  utilizing  the  two -wave length 
(4880R  and  10.  6  p)  propagation  facility  and  an  array  of  microthermal  sensors, 
the  latter  for  spatial  correlation  measurements  on  the  intermittent  turbulence 
field.  The  raw  data  are  magnetically  recorded  and  then  processed  on  a  PDP-11 
computer,  which  is  programmed  to  perform  a  variety  of  functions: 

Probability  Computations 

Probability  distribution  of  scintillation  and  single- 
or  double-probe  microthermal  data 

First  four  moments  and  central  moments 

Skewness,  kUrtosis 

Above  quantities  for  absolute  value  of  signal  (with 
corrected  mean) 

Above  quantities  for  square  of  signal  (or  other 
function  of  signal) 

Above  quantities  with  recursive,  first-order  digital 
(RC)  filter,  for  signal  or  function  of  signal 
(variable  averaging  before  computation  of 
probability  parameters) 
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Above  quantities  for  first  or  second  (finite -difference) 
derivatives  of  signal 

Variance  vs  averaging  time  of  signal  or  function  of 
signal,  and  mean  and  variance  of  variances 

Spectral  Computations 

FFT  for  spectra  from  1  Hz  to  500  Hz  with  1  Hz  resolu¬ 
tion  (averaged  over  many  individual  spectra  taken 
over  1  sec  each) 

Low-frequency  spectra  from  low-pass -filtered  signal 
(lowest  frequency  limited  by  length  of  data  record) 

Frequency-weighted  spectra 

Log  spectra. 

The  major  goals  of  this  investigation  will  be  accomplished  in  a 
follow-on  program.  In  this  report,  we  will  briefly  indicate  the  types  of 
data  and  analyses  which  are  pertinent,  and  the  theoretical  issues  involved. 

B.  Temperature  Statistics 

There  has  been  a  significant  amount  of  theoretical  work  in  the 
literature  on  the  statistical  properties  of  turbulent  fields,  including  hypo¬ 
theses  on  the  forms  of  the  distributions,  the  nature  of  the  moments,  and 
the  effects  of  averaging.  In  our  preliminary  investigations,  we  have  found 
better  empirical  agreement  with  certain  of  these  theoretical  treatments  than 
has  been  heretofore  reported.  The  details  are  described  below. 
Microthermal  Fluctuations  and  Scintillations  (Qualitative) 

Typical  single -probe  fluctuations  for  strong -turbulence  conditions 
and  a  low  (  <  1  m/sec),  unstable  wind  speed  are  shown  in  Fig.  14.  Although 
this  particular  record  does  not  represent  conditions  of  pronounced  inter- 
mittency,  the  sharply  asymmetrical  nature  of  the  fluctuations  is  apparent. 

In  Fig.  14b^,  a  ramp-like  structure  may  be  discerned  as  described  in  the 
literature.  Similar  fluctuations  for  a  steady,  6  m/sec  wind  are  shown  in 
Fig.  15.  nven  though  these  records  for  highly  developed  turbulence  are  not 
highly  intermittent,  the  very  "spiky"  nature  of  the  process  shows  that  the 
temperature  (and  hence  refractive  index)  fluctuations  are  highly  nongaussian. 
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Under  conditions  of  weak  or  less-developed  turbulence,  such  as 
would  be  observed  on  a  cloudy  day  or  at  night,  the  intermittency  is  much 
more  pronounced  and  may  involve  relatively  long  spatial  or  time  scales. 

^n  example  of  this  is  shown  in  Fig.  I6r  The  similarly  sporadic  nature 
of  the  first  derivative  (first  increment)  is  also  shown.  However,  it  may 
be  noticed  that  the  optical  and  infrared  scintillations  do  not  manifest  this 
intermittency,  since  they  represent  (weighted)  averages  of  the  turbulence 
field  over  the  path.  However,  in  the  more  interesting  cases,  the  inter¬ 
mittency  is  very  pronounced  on  the  scale  of  a  few  seconds,  8  and  the 
optical,  infrared  scintillations  also  manifest  a  degree  of  sporadicity. 

Double-probe  fluctuations  inherently  filter  out  the  low  frequency 

temperature  components  (corresponding  to  spatial  scales  greater  than 

the  probe  separation),  and  are  symmetric  (Fig.  17).  In  the  presence  cf 

substantial  intermittency,  this  random  process  can  be  described  by  a 

fast,  symmetrical,  spiky  process  of  zero  mean  multiplied  by  a  slow, 

1  8 

random  envelope  process.  '  This  will  be  discussed  further  below. 
Probability  Distributions  of  Microthermal  Fluctuations 

Very  little  can  be  said  concerning  the  probability  distributions  of 
the  raw,  single-  or  double-probe  microthermal  fluctuations  (Figs.  14-17), 
except  that  the  kurtosis  is  large.  The  distributions  for  single  and  double 
probes  are  shown  in  Fig.  18  a,  b  respectively;  the  single-probe  distribution 
is  asymmetrical  as  expected,  while  that  of  the  double  probe  is  basically 
symmetrical.  In  18b,  a  normal  distribution  with  the  same  variance  is 
shown;  the  actual  distribution  is  much  more  peaked  around  zero,  and  from 
the  value  of  the  kurtosis  (21.3),  the  probability  of  large  values  is  seen  to 
be  much  greiter  than  for  the  gaussian  case. 

Envelope  Probability  Distributions 


In  the  case  of  pronounced  intermittency,  the  probability  distribution 
of  the  envelope  function  is  of  interest.  1,8  In  order  to  develop  this  distri¬ 
bution,  the  absolute  value  of  the  two -probe  fluctuations  (  |  A  T  |  )  is  taken, 
and  then  the  signal  is  smoothed  over  an  averaging  time  of  a  few  seconds. 
This  results  in  the  spectral  components  of  the  multiplicative  envelope 


function  being  de  /eloped  as  additive  low  frequencies  rather  than  simply 
as  narrow  sidebands  about  the  high-frequency  fluctuations.  The  resulting 
distribution  is  approximately  indicative  of  a  bimodal  distribution,  35  so 
that  the  simplest  model  of  intermittency  would  be  described  by  two  discrete 
levels  of  turbulence  with  random  transitions  between  them.  As  the  averaging 
time  is  increased  to  many  decorrelation  times  for  the  envelope,  the  distri¬ 
bution  of  course  becomes  gaussian. 

"Good"  Turbulence  Variables 


The  temperature  field  is  not  a  "good"  characteristic  of  the  turbulence 
because  it  is  neither  stationary  nor  homogeneous.  However,  it  does  have 
stationary  first  increments,  8  so  that,  in  particular,  its  spatial  (or  tem¬ 
poral)  derivative  is  stationary.  This  will  also  be  true  of  two-probe  differ¬ 
ential  measurements,  as  long  as  the  probe  separation  is  not  too  large. 

These  derivatives  and  finite  differences  are  therefore  "good"  variables. 

Theoretical  treatments  have  predicted  that  non-negative,  "good" 
variables  will  exhibit  log  normal  behavior  if  averaged  over  a  spatial  ( jr 
corresponding  temporal)  scale  within  the  inertial  subrange.  A  related 
experiment  in  the  literature  on  the  squared  temperature  derivative  showed 
a  slight  skewness,  which  was  attributed  to  lack  of  local  isotropy.  38 

The  squared  first  derivative  is  shown  in  Fig.  19  for  various  averaging 
times.  The  probability  distribution  is  shown  in  Fig.  20,  and  is  seen  to  be 
very  similar  to  that  of  Fig.  18b,  with  a  kurtosis  of  25.  6.  However,  the 
distribution  of  the  absolute  value  of  the  derivative  vs  averaging  time  is 
shown  in  Fig.  21,  and  it  is  seen  that  (over  the  probability  range  covered) 
good  log  normality  is  exhibited  for  10  msec  averaging.  Since  the  wind  speed 
for  this  case  was  6  m/sec,  the  spatial  averaging  scale  was  6  cm--which  is 
comfortably  within  the  inertial  subrange.  Although  the  process  retains  its 
"log-normal-like"  nature  for  greater  and  less  averaging,  noticeable 
departures  from  strict  normality  are  observed. 


A  useful  way  t r  characterize  a  random  process  is  on  a  "beta 
diagram",  34  which  is  a  plot  of  f^2)  vs  (p,2/^3,,  whore  %  is  a  central 
moment  of  order  n.  This  plot  is  shown  vs  averaging  time  in  Fig.  22,  for 
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the  absolute  value  of  the  temperature  derivative.  The  solid  line  indicates 
the  locus  of  log  normal  distributions,  and  the  point  represents  a  gaussian 
distribution.  It  is  seen  that,  as  averaging  increases,  the  points  depart 
somewhat  from  the  log  normal  line  and  converge  to  the  gaussian  point  as 
would  be  expected.  The  departure  from  log  normality  represents  an  inter¬ 
esting  and  subtle  manifestation  of  the  outer  scale  of  turbulence. 

37 

The  theory  of  Gurvich  and  Yaglom  and  the  properties  of  log  normal 
distributions  indicate  power  law  behavior  of  the  moments  as  functions  of 
averaging  time  (t  ): 

-a 

(jl  ~  T  n  (28) 

n 


where  T  corresponds  to  spatial  scales  within  the  subrange.  A  plot  of  the 
first  four  moments  of  the  distribution  of  Figs.  21  and  22  vs  averaging  time 
is  shown  in  Fig.  23.  The  power  law  behavior  is  confirmed  for  scales  within 
the  subrange;  however,  the  values  of  the  exponents  do  not  agree  with  the 
theoretical  predictions. 

C.  Relationship  Between  Optical /infrared  and  Microthermal 
Quantities 

As  described  in  Refs.  1  and  8,  a  major  goal  of  this  effort  is  to  relate 

data  spread  and  confidence  statements  for  short-term-average  optical/infra- 

red  measurements  to  corresponding  microthermal  quantities.  To  the  extent 

that  pronounced  intermittent  enters,  the  interesting  averaging  time  for  the 
2 

latter  (C  )  is  on  the  order  of  a  few  seconds,  so  that  fast  fluctuations  are 

n  2 
removed  blithe  intermittent  envelope  is  not.  As  explained  above,  Hie 

distribution  may  be  quasi-bimodal,  and  any  intermittency  effects  showing 

up  in  the  scintillations  must  result  in  a  distortion  of  the  log  normality  of 

irradiance. 

Work  is  continuing  on  the  necessary  modelling  and  theoretical  rela¬ 
tionships.  An  empirical  determination  of  log  amplitude  (scintillation) 

2  2 

variances  (<r_  )  vs  spread  in  microthermal  C  measurements  is  shown  in 
E  2  n 

Fig.  24.  The  individual  values  were  taken  for  a  period  of  10  seconds,  and 
the  spread  and  mean  of  a  number  of  such  values  is  shown  for  e  ich  of  six  runs. 


Since  the  scintillations  were  measured  at  10.  6  p.,  saturation  of  scintilla¬ 
tions  was  not  a  factor;  the  line  indicates  agreement  between  experimental 
and  theoretical  variances,  where  the  latter  is  given  by^ 


<rT2=  0.  124  C 


W1/6 

n 


(29) 


Similar- -but  much  smaller --error  bars  could  also  be  constructed  from  a 

2 

sequence  of  measurements  of  a-  for  each  run. 

E 

It  is  apparent  from  Eq.  (29)  that  an  analytical  consideration  of  the 
spread  obtained  in  measurements  of  optical/infrared  variances  will  involve 
3patial  correlations  of  short-term  microthermal  variances  (or  C  2 ): 

Var  <r  2  --  /c  2  C  2 
E  \  "1  "2 

For  this  reason,  arrays  of  portable  probe  systems  are  being  employed. 

D.  Realizations  of  Random  Temperature  Field  and  Monte  Carlo 
Simulations 

In  addition  to  analytical  and  empirical  approaches,  we  intend  to 
generate  realizations  or  ensemble  members  of  the  random  temperature 
field,  and  to  employ  Monte  Carlo  simulations  of  the  propagation  problem. 

E.  Physical  Basis  of  Intermittency 

The  physical  basis  of  turbulence  intermittency  was  recently 

3  9 

reviewed  by  Mollo-Christensen,  and  the  reader  is  referred  to  that 
paper.  Turbulence,  including  the  geophysical  case,  is  described  as 
occurring  sporadically  or  in  "bursts",  which  are  episodes  of  generation 
involving  energy  in  a  nonlinear  interaction  through  a  hierarchy  of  (coupled) 
scales.  These  different  scales  of  instability  synergize - -i.  e.  ,  they  inter¬ 
act  to  mutually  enhance  growth  by  extracting  energy  from  the  large  scales. 
The  higher  the  Reynolds  number,  the  shorter  and  more  intensive  the  bursts. 
The  physical  arguments  suggest  fruitful  areas  of  microthermal  experimenta¬ 
tion,  such  as  the  simultaneous  joint  measure  of  several  scales. 
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Of  great  importance  in  the  present  context  is  the  universality  of 
intermittency  in  turbulent  flows;  the  effect  of  topography  and  altitude;  and 
the  relationship  between  fundamental  intermittency  as  described  here,  and 
phenomena  such  as  plumes. 

F.  Related  Theoretical  Issues 

The  basic  model  of  the  intermittent  thermal  field  in  terms  of  a 

multiplicative,  low-frequency  function  describes  a  random  process  which 

1  8 

has  certain  unusual  properties.  In  particular,  ’  the  power  spectrum  and 
autocorrelation  of  the  process  do  not  exhibit  the  sporadicity,  and  e.  go 
averaging  time  analyses  based  on  these  quantities  are  invalidated.  This 
raises  the  issue  of  whether  or  not  such  a  process  requires  a  more  general 
approach  from  stochastic  theory,  and  several  aspects  have  been  examined 
as  follov  s: 

Orthogonal  Increments 

The  fact  that  the  low-frequency,  random,  multiplicative  envelope 

establishes  narrow  sidebands  around  each  (continuum)  component  of  the 

basic  microthermal  spectrum  raises  the  question  of  suitability  of  the 

F ourier -Stieljes  representation,  and  the  orthogonality  of  the  spectral 
36  , 

increments.  We  have  shown  that,  providing  the  mu.  tiplicative  function 
is  stationary  in  the  wide  sense,  the  overall  process  has  orthogonal  incre¬ 
ments.  This  conclusion  will  be  subjected  to  an  empirical  test  on  the  data. 
Generalized  Stationarity 

In  the  above  model  of  intermittency,  there  exist  no  additive,  low- 
frequency  components  in  the  power  spectrum  of  the  actual  microthermal 
process.  However,  such  components  are  often  physically  present--e.  g.  , 
in  plumes,  which  are  warmer,  more  turbulent  regions  relative  to  their 
surroundings.  Such  low-frequency  components  can  cause  serious  problems 
in  stochastic  theory,  and  in  fact  first  increments  such  as  differential 
probe  measurements  are  often  utilized  to  avoid  such  stationarity  questions. 

The  behavior  of  a  random  process  at  very  low  frequencies  (the 
"infrared  catastrophe")  is  treated  by  Mandelbrot,  40  in  which  he  generalizes 
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Jie  concept  of  stationarity  to  conditional  statements,  and  treats  "sporadic 
processes".  However,  since  scintillation  effects  relate  to  turbulence 
components  generally  in  the  inertial  subrange  (hence  the  applicability  of 
differential  rather  than  single  probe  measurements),  we  believe  that 
these  low  frequency  components  are  irrelevant  to  our  problem.  This  may 

not  be  entirely  the  case  for  optical  infrared  phase  effects,  which  involve 
larger  atmospheric  scales. 

Order  and  Dimension  of  a  Random  Proems36 

The  "order"  of  a  random  process  refers  to  the  order  of  derivative 

or  finite  increment  necessary  for  stationarity.  As  discussed  above,  the 

basic  nature  of  the  thermal  field  necessitates  the  use  of  first  increments. 

However,  it  is  not  apparent  that  intermittency  relates  to  higher  "orders" 
as  here  defined. 

The  "dimension"  of  a  process  refers  the  dependence  of  e.  g.  struc¬ 
ture  functions  on  independent  variables  which  are  tensors.  Again,  it  is 

not  apparent  that  this  relates  to  intermittency,  although  it  does  relate  to 
anisotropy. 

G.  The  Scintillation  Process:  Conditional  Statistics,  Level 
Crossings,  Deep  Fades,  Markov  Properties 
We  are  interested  in  the  conditional  probatility  of  the  instantaneous 
irrad lance  at  time  t.,  given  the  value  at  t, ,  and  related  issues.  Unfortunately 
we  have  empirically  found  the  process  to  be  non-Markovian.  To  the  extent 
that  the  log  amplitude  is  gaussian,  however,  the  problem  is  readily  tractable. 
Also,  in  order  to  attack  the  problem  of  deep  fades  in  an  intermittent  type  of 
process,  use  may  be  made  of  the  theory  of  extreme  values.41 
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1.  Asymptotic  Log  Curves  for  Total  Angular  Beamspread  vs  (b/p  ). 

"  o 

a.  fixed  and  b  varied. 

b.  b  fixed  and  u  varied. 

o 

2.  The  Function  f(b/p  )  Inferred  from  the  Numerical  Results  of  Refs.  16 
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a.  vs  Log  b  for  fixed  p  ,  ,  and  c r  2 

ro  v  K  ’ 


b.  vs  Log  p 
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6.  Power  Spectrum  (S)  and  Frequency-Weighted  Power  Spectrum  of 

Transmitter  Tracking  (Wander)  Signal. 

7.  Cumulative  Probability  Distribution  of  Vertical  Component  of  Trans¬ 

mitter  Tracking  (Wander)  Signal. 

8.  Target  Irradiance  at  Beim  Axis  with  1  Second  Averaging 

a.  Without  Wander -tracking 

b.  With  Wander -tracking  . 


9.  Measured  Mean  Target  Irradiance  vs  (b/pj,  With  and  Without  Tracking. 

10*  Rati°  ^Measured  Mean  Target  Irradiances  With  vs  Without  Tracking, 
o 

11.  Power  Spectra  (S)  of  Log  Amplitude  Irradiance  (Fading)  With  and  Without 

Tracking  (b/p  =113). 

12.  Frequency  Weighted  Power  Spectra  of  Log  Amplitude  Irradiance  With 

and  Without  Tracking  (b/p  =  68). 

13.  Cumulative  Probability  Distribution  of  Log  Amplitude  Irradiance  With 

and  Without  Tracking  (b/p  =  96). 

14'  Sing^;^roobe  Microthermal  Fluctuations  for  Strong  Turbulence  and  Low 
Wind  Speed  (  <  1  m/sec  ) 

a.  10  sec/div 

b.  1  sec/div. 

15.  Same  as  Fig.  14,  for  6  m/sec  Wind  Speed. 

16.  Microthermal  and  Optical/infrared  Fluctuations  at  Night  (Weak 

Turbulence) 

a.  Single -Probe  Fluctuations 

b.  Square  of  Derivative  of  (a) 

c.  Log  Scintillations  at  4880  R 

d.  Log  Scintillations  at  10.6  p. 

17.  Double-Probe  Microthermal  Fluctuations  for  Strong  Turbulence  and  l  ow 

Wind  Speed  (  <  1  m/sec), 

18.  Probability  Density  F unction  for  Microthermal  Fluctuations 

a.  Single -Probe  Fluctuations  T 

b.  Double -Probe  Fluctuations  A  T.  The  Gaussian  Distribution 

With  the  Same  Variance  is  Also  Shown. 
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19.  Squared  Temperature -Derivative 


a.  1  sec  Averaging  Time 

b.  0.  1  sec  Averaging  Time 

c.  0.  01  sec  Averaging  Time. 

20.  Probability  Density  Function  for  Temperature  Derivative.  The 
Gaussian  Distribution  With  the  Same  Variance  is  Also  Shown. 


21.  Cumulative  Probibility  Distribution  for  Log  of  Absolute  Value  of 

Temperature  Derivative,  for  Different  Averaging  Times. 

22.  (3-Diagram  for  Absolute  Value  of  Temperature  Derivative,  for 

Different  Averaging  Times.  The  Unaveraged,  Linear  Value  is 
Also  Shown. 

23.  First  Four  Moments  of  Absolute  Value  of  Temperature  Derivative, 

vs  Averaging  Time. 

24.  Experimental  Log  Amplitude  Variance  of  10.  6  ^  Scintillations  vs 

Theoretical  (Rytov)  Value. 

Each  Theoretical  Value  is  Based  on  an  Individual  C  2 
Measurement  Taken  Over  a  10  sec  Averaging  Time.**1 
The  Standard  Deviation  and  Mean  of  a  Number  of  Such 
Measurements  is  Shown  for  Each  of  Six  Runs. 

Runs  l,  3,  5:  Low  Wind  Speed 
Run  6:  Essentially  Zero  Wind 
Runs  2,4:  Wind  Speed  >  4  m/sec. 
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